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Abstract 



We study models of chiral interacting fermions by means of conformal 

and Bethe-Ansatz techniques, and determine their thermodynamic properties 

and asymptotic correlation functions. We identify a class of fixed points 

characterizing the infrared behavior of the models. They display chirally 

stabilized non Fermi-Liquid behavior characterized by universal exponents. A 

realization of these fluids may be found in the edge states of bilayered QHE 

systems. We calculate the universal temperature dependence of the staggered 

conductance in these systems. 
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Models displaying Non-Fermi-Liquid (NFL) behavior have been intensely studied follow- 
ing the suggestion that the normal phase of the cuprate superconductors cannot be described 
in terms of the conventional Landau theory The study of such models in d > 2 has made 
slow progress. In one dimension, however, a class of models was discovered, the Luttinger 
Liquids, exhibiting NFL behavior characterized in the infrared (IR) limit by single particle 
correlation functions having only cuts with non-universal exponents. 

In this letter we identify another large class of IR fixed points, available to interacting 
fermions with unequal number of left and right moving degrees of freedom. The new fixed 
points control the spin and (or) flavor sector of the theory and are universal in character, 
independent of the strength of the interaction. We shall refer to these fixed points as chirally 
stabilized or coset fluids. 

Several systems may be described by these chiral fluids. They must break T-invariance, 
either by the presence of a magnetic field or by the interactions. Below we shall discuss in 
some detail one example, the bilayered Quantum Hall system, and present predictions for 
the temperature dependence of its staggered conductance. 

The models we study tend in the ultraviolet (UV) to the free hamiltonian, 

H = -Wpfdx (j2 f r=l^R, a ,r( X ) d ^R^A X ) ~ ££l ^L^l^Wx^ L,a,l{ X )) • The fields ^*R,a,r{ X ) 

(ipLai( x ))j whh a = ±1 the spin index, and r = 1.../r, (I = the right (left) flavor 

index, create right-(left-) moving particles with a linearized dispersion e = ±vp{k — hp)- 
These particles are conventionally considered in the Fock basis generated by the Fourier 
modes of i>*R ar {x) and i'l^x). As such H Q is an example of a (non-interacting) Fermi 
liquid. 

The kinematics of one dimension permits many other bases to describe particles with 
linearize dispersion, since a linear combination of any number of left (right)-movers is again 
a left (right)-mover. This also manifests itself in the operator language, in the representation 
of the field via abelian (or non-abelian) bosonization allowing a separation of the charge, 
flavor and spin components. The abelian basis U(1)^ L (g> U{ly R is convenient when the 
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interaction terms involve the charge densities, = J2i Pl,i(%) = J^a,^! aii^^L^iix) 

(similarly for p R ). In this basis the fermionic fields take the form, ipL,a,i — e l ^ L > a ' 1 ^ with 
<t>L,a,i{.x) a left moving bosonic field. The non-abelian basis (S77(2)/ L x SU(fL)2 x U{1)) ® 
(SU(2)f R x SU(fn)2 x U(l)) is convenient when the interactions are mediated by the spin- 
densities Si,S R , where S l L (x) = EiS l L /x) = E^i^l^iK^lai (similarly for S R ), which 
close on Kac-Moody algebras with central charges /x and respectively. Hq is then 
represented as a sum of theories @FJ] corresponding to the charge- spin- and flavor- compo- 
nents Q, and the operators of the free fermi theory are expressed in terms of the operators 
in the component theories, ipR, a ,r(x) = 9r{ x ) h r R (x) e l $ R ^ where g a and h r are fields in 
the spin and flavor sectors, respectively, transforming in the fundamental representation of 
SU(2) and SU(f R ), and <p R is the charge field. Each component of Hq has its own charac- 
teristics. In particular, the spin component, SU{2)f L ® SU{2)f R , has the Virasoro charge 
c L + c R = 3f L /(f L + 2) + 3f R /(f R + 2) and the Kac-Moody charge k L + k R = f L + f R . 
Other bases are available: applying the Bethe Ansatz technique to the model allows the 
construction of an arbitrary number of bases corresponding to a choice of a scattering ma- 
trix between the left- (right-) movers Sll(S rr ) Still the model describes a FL since the 
various components can be put together to form a fermion. 

When interactions are added the situation can change. To lead to a new behavior in the 
IR the interaction needs to flow to some new fixed point, preferably at intermediate coupling 
since strong coupling fixed points tend to open a gap. The different components making up 
the electron are then sufficiently modified, and one may find that in the infrared the electron 
can no longer be reconstituted. In other words, the natural degrees of freedom will no longer 
be fermionic and one will observe a NFL-behavior. This behavior is manifested in the 
structure of the single particle correlation functions at large distances, less so in the higher 
correlation functions and the low temperature thermodynamics. Thus the specific heat of 
any model characterized in the IR by some conformal fixed point hamiltonian will be will 
be linear, cy = ^{cl + cr)T, and the magnetic susceptibility is a constant x — (kz + k R )u Q , 
(i/ = 1/ttvf the density of states per unit length) whether the fixed point describes a FL or 
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not. 

A familiar example is the Luttinger Model, obtained by modifying the charge sec- 
tor. The hamiltonian is (here we choose fa = /&), H = H + g 2c J dxpL(x)pn(x) + 
<74 C / dx (pl{x)pl{x) + pr(x)pr(x)) . Renormalization group calculations || indicate that the 
model is conformally invariant (in the infinite cut-off limit.) As is well known, upon express- 
ing the hamiltonian in terms of bosonic fields it becomes quadratic and can be solved by a 
Bogoliubov rotation. The g± term modifies the spin and charge velocities (without destroy- 
ing the FL property at the fermi surface [0]), while the effect of the #2 term is to modify the 
exponents of the fermionic correlation functions destroying the pole structure characteristic 
of a FL. Many models - the Luttinger liquids ||- exhibit this low-energy behavior: they flow 
in the infrared to a c = 1 bosonic model characterized by continuous parameters with no 
fermionic interpretation. 

We proceed now to study models flowing in the IR to fixed points describing the class of 
universal non Fermi-liquids. The models are obtained by adding spin-exchange interactions 
to Hq, 

H = H + J2jdx (gsfS^Sl^x) . (1) 

r,l 

with (g s ) a matrix of couplings. The models are isotropic when the matrix is diagonal. 

Standard calculations show that the perturbation destabilizes the weak coupling fixed 
point. To identify where the models flow to one may perform a strong coupling or a large- 
(/x+/r) expansion. Instead we shall apply techniques from conformal field theory and 
thermodynamic Bethe Ansatz (TBA) and conclude that the models flow to fixed points, the 
chirally stabilized fluids, given by a particular class of conformal field theories - the WZW 
coset models (see below). We shall show that they describe NFL behavior by studying the 
correlation functions and the S-matrices. 

Here we chose an interaction term modifying the spin sector only. Other interactions 
may be added. Terms of the form g c J dxpn(x)pL(x) + / dx (gf) ab F^ a (x)F^ b (x) (F^ a = 
^Rar^rr'^R^y is the flavor density) will modify the charge and flavor sectors: the g c coupling 
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will turn the isolated fixed point into a line of NFL- fixed points, while the gf term will drive 
the model to a non-trivial fixed point in the flavor sector. One may further add terms of 
the form / (S l L (x) S l L (x) + S R {x)S R (x)) or / (f£(x)F£(x) + F&(x)F&(x)) which will modify 
the spin and flavor velocities, respectively, or also terms that break global spin (or flavor) 
invariance. These issues will be discussed elsewhere. 

We begin by discussing the flavor-isotropic model, characterized by one coupling g s . The 
interaction breaks the spin symmetry to a global SU (2). Therefore it is no longer conformally 
invariant, its IR properties will depend on Jr-Jl- For /r=/l, the model is chirally invariant 
and flows to a strong coupling fixed point generating a mass gap. For J'r > Jl, on the other 
hand, the model is gapless and flows to a non-trivial fixed point which we proceed to identify. 
This is possible since in this case the model is chiral in a strong sense: the conformal central 
charges as well as the Kac- Moody central charges on left and right are different. The charge 
differences cr-cl and must, however, be preserved under the flow! 0. This places a 

strong constraint on the IR fixed point - there is a unique fixed point theory of lowest central 
charge satisfying these two conditions. Thus, 

SU(2) fL ® SU(2) fR — SU{2 \^y«- fL ® SU(2) fR „ fL . (2) 

Here the Kac-Moody central charge f R — fi is matched by postulating that the right 
degrees of freedom are those of the chiral WZW model SU(2) f R -f L while the left degrees of 
freedom are singlet in spin and are described by chiral coset CFT |10| . 



The specific heat and the magnetic susceptibility can be immediately determined from 
the IR central charges of the Virasoro and Kac-Moody algebras in the IR theory. Hence 
the specific heat will be linear in the temperature in the UV and in the IR limits (with 
corrections |PH|), and will undergo the flow: 



c v --(h + f R )T -> c v --{f L + f R + fR _ h + 2 -—)T 

where we also included charge and flavor contributions. The flow in the suceptibility 
will be, x uv = Un + fi)vo —> X ir = {fa ~ h)vo, leading to a Wilson ratio R w = 



(/l + fa + fx-fL+2 ~ fi+2) Kf R ~ f L ^' ^° determine the scales where the crossover in 
behavior occurs one needs to construct the complete theory, this is done below through a 
Bethe-Ansatz. 

We discuss now the operators around the IR fixed point. A SU(2) k theory contains 
primary fields (fP m ^ m {x,t) transforming under a particular left and right representation of 
the symmetry There is a finite number of operators allowed: < j,j < k/2, and their 
dimension is h — ~ ? ^ 2 1 ^ . For the coset theory there is a single primary ^-9, for each choice 
< j < /l/2, < f < (f R - f L )/2 and < j" < f R /2. The dimension of the primary 
is the difference of the dimensions of the group primaries, up to an integer. We can thus 
match the physical fields with the operator basis around the fixed point and read off the IR 
behaviour of the correlation functions, 

<^,a,,(M)^w(0,0) > ^5 aa '5 ll \x-v F t)-^\x + v F t)- 5 - 

< Si(x, t)S{{0, 0) > -> 5 ij (x + v F t)~ 2 (x 2 - v 2 F t 2 y 4/(fR+2) 

< S l R (x,t) 4(0,0) > -> 5 lj {x-v F t)- 2 

with 5 L = 3/2(/ fl +2) and 8 R = 3f L /2(f R -f L +2)(f R -f L ). We observe that the FL structure 
is destroyed: the momentum distributions for small momenta are n a (k) ~ \k - k F \ 2Sa . In 
these expressions the left and right components move with velocity v F . The inclusion of the 
term / (J l L J l L + J R J R ) would modify this. Also, the charge and flavor correlation function 
will be non-trivial upon inclusion of the terms mentioned earlier. 

We consider now the case with flavor anisotropy. Begin by studying the various limits of 
extreme anisotropy, which can be modeled as a sequence of flows, each of the type described 
above. Consider for example a coupling g rl = g 1 for r < f R i and g rl = g-i for r > f R1 . This 
breaks the SU(f R ) x U(l) right flavor and charge symmetry down to SU(f R1 ) x f/(l) x 
SU(f R2 ) x U(l) with f R1 + f R2 = f R . Clearly we want to bosonize the two groups of right 
fermions separately, introducing spin densities S R1 and S R2 generating SU{2) Kac-Moody 
algebras of level f R \ and f R 2, so that the interaction will again involve only the spin sector 
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of the theory. 

In the limit g\ ^> g 2 , we can regard the g\ interaction as generating precisely the flow 
described above, approaching arbitrarily closely to the IR fixed point described above. We 
can then identify the g 2 interaction as a specific perturbation of this IR fixed point using 
our earlier results. If it is still marginally relevant, this will produce a flow to a final IR fixed 
point. Of course this analysis would be reversed for g 2 <C g%. For the intermediate regime, 
we can make a guess as to the likely behavior by appealing to the c-theorem: non-trivial 
flows in 1 + 1 dimensions always decrease c | 12|| . If we compare the two final IR fixed points 



reached by the two limits of extreme anisotropy and find that one has higher c, it is likely 
that any finite anisotropy will cause the flow to continue to the other IR fixed point. 

There are several patterns which can arise in our example. If fm > fi > 
f R2 , the gx > g 2 limit will start with the flow SU{2) fL ® SU(2) fm x SU(2) fm -> 

f su(2) f — fm ~ L ® SU(2)f R1 _f L x SU(2)f R2 . The remaining interaction is irrelevant at this 

fixed point. The g 2 ^> g\ limit will follow a different sequence: SU(2)f L £g> SU{2)f R2 x 
SU(2) fm - SU{2) fL _ fm ® Sm ^^-^ x SU(2) fm - ^)/,-/^x5^ /m+/m - / , ^ 
su(2) fR2 ^xsu(2) fL _ fR2 x SU(2) fR1+fa2 _ fL . If fi > f m > f R2 , the two limiting sequences are 
both of the latter form - precisely this if g 2 ^> g%, and with fm and f R2 interchanged for 
9\ 92- One can check that in either case, if fm > f R2 , the result of the g 2 3> g\ sequence 
always has lower c than the result of the g\ ^> g 2 sequence, making it the IR fixed point for 
generic anisotropy. The correlation functions for this case can be found by the same means 
and will be given in a subsequent work. 

We can actually follow the flow at any scale by solving the model exactly. The model 
is closely related to the multichannel Kondo model and exhibits dynamical fusion |13| (a 
different approach was given in ||14|) allowing a solution by a method very similar to the 



one used to solve the anisotropic multichannel Kondo model We find that the model 
generates scales m l L , m r R : I < fi, r < fn, parametrizing the patterns of flavor symmetry 
breaking, and setting the excitation energies and momenta. The free energy is given by, 
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TL r°° ( h h \ 

F(T, h) = -— m r R e-t In(l + -)) + £ m^e^ ln(l + r^C, ^ )) J , 

where the functions {r/(£, ^)} are the solution of the following system of coupled integral 
equations (TBA-equations): 

\nr] n = -2^fe«-2^e- 5 + Gln(l + r ?n+1 ) + Gln(l + 7fc_i), n=l,...,oo, V7q = 0, 

with boundary condition: ln?y n — > 2rifih/T. The integral operator G is defined by the 
kernel l/(27r cosh(£' — £)). In the isotropic case, = mLdf^j, m r R = m#5/ flir . The 
forms of the driving terms m^e^ and rriRe~^ are characteristic of massless left and right 
moving excitations, and when both occur at the same level a driving term mcosh£ results 
indicating a mass gap. To be more explicit, consider the case of two-channels of right movers 
and one-channel of left movers. The two couplings g±, g2, are obtained by diagonalizing 
the matrix of couplings g r s l . Choose g% < gi, 4> = g\J ' g%- The physical scales then are, 
= 2_D R cos( Z 7p)e~»2 ; m 2 R = D R e~~^ and m\ = 2D L e~~^ , with D Ll D R the densities of 



left and right movers [15 



From the TBA-equations both the IR and the UV limits can be read off using TBA- 



rules ||16|| : the IR limit of the left movers is obtained from the equations by considering the 
right-mass as infinitely heavy and truncating the equations at the level it was inserted. An 
analogous rule holds for the right movers. The UV limit is obtained, on the other hand, 
by considering the masses as vanishing. Applying these rules we deduce the IR limit, and 
find accord with the conformal considerations. The solution of the equations provides in 
addition the full interpolation between the UV and IR limits. 

To observe chiral NFL behavior in an experimental system one may study the edge states 
in a bilayer of a quantum Hall liquid. In the higher hierarchy states (e.g. for filling factor 
v = n/(np + 1) with p a negative even number) one edge state mode moves in a direction 
opposite to the rest []T7| providing a chiral imbalance. The channel degrees of freedom play 
the role of the flavor and the bilayer provides a 'spin' degree of freedom. The spin-exchange 
interaction (eq.(P) will be induced through virtual hoppings between the layers. We assume 
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that the couplings have been adjusted to have spin SU(2) symmetry. The system will then 
exhibit a universal behavior in the staggered conductance G s measuring the response to 
electric fields that are oppositely oriented in the layers. The staggered current is given by 
J s = v f(S\ — Si), - it is the spin current of the model - and the staggered current-current 
correlation function can be expressed in terms of the spin density correlations calculated 
above. We conclude (setting fi = 1, fn = n — l) that Gs ~ 2 1 8/(™+ 1 )- 1 ) an( j expect for 



n = 7 a staggered metal-insulator transition induced by correlations fl8 
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